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(Dated: April 16, 2007) 

We describe a search for the standard model Higgs boson with a mass of 105 GeV/c^ to 145 GeV/c^ 
in data corresponding to an integrated luminosity of approximately 450 pb~^ collected with the DO 
detector at the Fermilab Tevatron pp collider at a center-of-mass energy of 1.96 TeV. The Higgs 
boson is required to be produced in association with a Z boson, and the Z boson is required to 
decay to either electrons or muons with the Higgs boson decaying to a bb pair. The data are well 
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described by the expected background, leading to 95% confidence level cross section upper limits 
a{pp ZH) X B{H -> bb) in the range of 3.1 pb to 4.4 pb. 

PACS numbers: 13.85.Ni, 13.85.Qk, 13.85Rm 



Over the past two decades, increasingly precise exper- 
imental results have repeatedly validated the standard 
model (SM) and the relationship between gauge invari- 
ance and the embedded coupling strengths. For massive 
W and Z bosons, gauge invariance of the Lagrangian is 
preserved through the Higgs mechanism, but the Higgs 
boson (H) has yet to be observed. The current lower 
bound on the mass of the Higgs boson from direct ex- 
perimental searches is Mh = 114.4 GeV/c^ at the 95% 
confidence level [l[. Searches for pp WH —^ e{fi)i/bb, 
pp WH WWW* a,nd pp ^ ZH ^ wbb have 
been recently reported [3, i, ij. The CDF collabora- 
tion has reported results in the pp — > WH — > and 
pp ZH £^£^bb {£ = e,u) channels with signifi- 
cantly smaller data sets 0, H, 0| ■ This Letter provides 
the first results from the DO experiment of searches for 
a Higgs boson produced in association with a Z boson, 
which then decays either to an electron pair or to a muon 
pair. The Higgs is assumed to decay to a bb pair with a 
branching fraction given by the SM. The Z{—^ i^£~)H 
channels reported in this letter comprise major compo- 
nents of the search for a Higgs boson at the Tevatron 
collider. 

Z bosons are reconstructed and identified through 
pairs of isolated, electrons or muons with large momen- 
tum components transverse to the beam direction {pr) 
having invariant mass consistent with that of the Z bo- 
son. Events are required to have exactly two jets iden- 
tified as arising from b quarks (6 jets). The resulting 
data are examined for the presence of a {H bb) sig- 
nal in the ^-tagged dijet mass distribution. An efficient 
6-identification algorithm with low misidentification rate 
and good dijet mass resolution are essential to enhance 
signal relative to the backgrounds. The analysis of the di- 
electron [§] (dimuon channel is based on 450±27 pb~^ 
(370 ± 23 pb~^) of data recorded by the DO experiment 
between 2002 and 2004. 

The DO detector [13, U] has a central-tracking sys- 
tem consisting of a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT), both located within 
a w 2 T superconducting solenoidal magnet, with de- 
signs optimized for tracking and vertcxing covering 
pseudorapiditics |r/| < 3 and |r/| < 2.5, respectively 
[i] = — ln[tan(0/2)], with 9 the polar angle relative to 
the direction of the proton beam). Central and for- 
ward preshower detectors are positioned just outside of 
the superconducting coil. A liquid-argon and uranium 
calorimeter has a central section (CC) covering pseudo- 
rapidities up to I77I « 1.1 and two end calorimeters (EC) 
that extend coverage to \r]\ « 4.2, with all three housed 
in separate cryostats |ll| . An outer muon system, cover- 



ing 1 77 1 < 2, consists of a layer of tracking detectors and 
scintillation trigger counters in front of 1.8 T toroids, 
followed by two similar layers behind the toroids [l^ . 
Luminosity is measured using plastic scintillator arrays 
placed in front of the EC cryostats [13] ■ 

The primary background to the Higgs signal is the as- 
sociated production of a Z boson with jets arising from 
gluon radiation, among which Z+bb production is an irre- 
ducible background. The other background sources are 
tt production, diboson (ZZ and WZ) production, and 
events from multijet production that are misidentified 
as containing Z bosons. The backgrounds are grouped 
into two categories with the first category, called physics 
backgrounds, containing events with Z ot W bosons aris- 
ing from SM processes: inclusive Z + bb production, in- 
clusive Z + jj production in which j is a jet without 
b flavor, tt, ZZ, and WZ events. This background is 
estimated from simulation as described below. The sec- 
ond category, called instrumental background, contains 
those events from multijet production that have two jets 
misidentified as isolated electrons or muons which ap- 
pear to arise from the Z boson decay. This background 
is modeled using control data samples and the procedure 
described below. 

Physics backgrounds are simulated using the leading 
order ALPGEN 14 1 and pythia [l^ event generators, 
with the leading order CTEQ5l [l^ used as parton distri- 
bution functions. The decay and fra gme ntation of heavy 
flavor hadrons is done via evtgen [17|. The simulated 
events are passed through a detailed DO detector simula- 
tion program based on GEANT 18| and are reconstructed 
using the same software program used to reconstruct the 
collider data. The ZH signal, for a range of Higgs masses, 
is also simulated using pythia with the same processing 
as applied to data. Determination of the instrumental 
background and the normalization of the physics back- 
grounds are discussed below. 

Candidate Z ^ ee events are selected using a combi- 
nation of single-electron triggers. Accepted events must 
have two isolated electromagnetic (EM) clusters recon- 
structed offline in the calorimeter. Isolation is defined as 



calorimeter energy within AR < 0.4 of the electron direc- 
tion and E^^^j' is the energy in the electromagnetic por- 
tion of the calorimeter within AR < 0.2 of the electron 
direction. Candidate electrons must satisfy / < 0.15. 
Each EM cluster must have px > 20 GeV/c and either 
|77dct| < 1-1 or 1.5 < |77dot| < 2.5, where rjdct is the pseu- 
dorapidity measured relative to the center of the detector, 
with at least one cluster satisfying |77dct| < 1.1. In ad- 
dition, the lateral and longitudinal shower shape of the 



(0.2). ,p(0.2) 



^j^j in which E^^^^^ is the total 
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energy cluster must be consistent with that expected of 
electrons. At least one of the two EM clusters is also 
required to have a reconstructed track matching the po- 
sition of the EM cluster energy. Events with a dielectron 
mass of 75 < M^e < 105 GeV/c^ form the Z boson can- 
didate sample in the dielectron channel. 

Candidate Z — > events are selected using a set 

of single-muon triggers. Accepted events must have two 
isolated muons reconstructed offline. The muons must 
have opposite charge, pr > 15 GeV/c, and \r]\ < 2.0 
with muon trajectories matched to tracks in the central 
tracking system (i.e., the SMT and the CFT), where the 
central track must contain at least one SMT measure- 
ment. In addition, the central tracks are required to have 
a distance of closest approach to the interaction vertex 
in the transverse plane smaller than 0.25 cm. Muon iso- 
lation is based on the sum of the energy measured in the 
calorimeter around the muon candidate and the sum of 
the pT of tracks within AR = ^(A0)2 + {A-qY = 0.5 of 
the muon candidate normalized by the muon momentum. 
The distribution of this variable in background multi- 
hadron events is converted to a probability distribution 
such that a low probability corresponds to an isolated 
muon. The product of the probabilities for both muons 
in an event is computed, and the event is retained if the 
product is less than 0.02. Accepted Z boson candidates 
must have the opening angle of the dimuon system in 
the transverse plane (azimuth) of Aip > 0.4, and invari- 
ant mass 65 GcV/c^ < M^^ < 115 GeV/c^ This mass 
range differs from that of the dielectrons because of the 
difference in resolutions of electron energies and muon 
momenta. 

After selecting the Z candidate events, we define a 
Z+dijet sample which, in addition to satisfying the Z 
candidate selection requirements, has at least two jets 
in each event. Jets are reconstructed from energy in 
calorimeter towers using the Run II cone algorithm with 
AR = 0.5 [l^ with towers defined as non-overlapping, 
adjacent regions of the calorimeter At^ x Acf) = 0.1 x 0.1 
in size. The transverse momentum of each jet is cor- 
rected for multiple pp interactions, calorimeter noise, out- 
of-cone showering in the calorimeter, and energy response 
of the calorimeter as determined from the transverse mo- 
mentum imbalance in photon-|-jct events 2^, 21 1. Only 
jets that pass standard quality requirements and satisfy 
Pt > 20 GeV/c and \t]\ < 2.5 are used in this analy- 
sis. The quality requirements are based on the pattern 
of energy deposition within a jet and consistency with 
the energy deposition measured by the trigger system. 

For the Z ee channel, the normalizations of the 
smaller tt, WZ and ZZ backgrounds are computed using 
simulated events and next-to-leading-order (NLO) cross 
sections. Trigger efficiency, electron identification (ID) 
efficiency and resolution correction factors are derived 
from comparisons of data control samples and simulated 
events. The background contributions from Z+bj 



and Z + bb processes are normalized to the observed 
Z-fdijet data yield reduced by the expected contribu- 
tions from the smaller physics and instrumental back- 
grounds. The relative fractions of the Z + jj, Z + bj and 
Z+bb backgrounds in the Z-|-dijet sample are determined 
from the acceptance and selection efficiencies multiplied 
by the ratios of the NLO cross sections for these processes 
computed using the mcfm [i^] program and the next-to- 
leading order CTEQ6m parton distribution functions. 
For the Z — > fJ.'^ channel, all physics backgrounds are 
determined using simulated events with NLO cross sec- 
tions applied. Trigger efficiency, muon ID efficiency, and 
resolution correction factors are derived from comparison 
of data control samples and the simulated events. 

Instrumental backgrounds in both channels are deter- 
mined by fitting the dilepton invariant mass distributions 
to a sum of non-Z and Z boson contributions. The Z 
boson lineshape is modeled using a Breit-Wigner distri- 
bution convoluted with a Gaussian representing detector 
resolution. The non-Z background, consisting of a sum 
of events from Drell-Yan production and instrumental 
background, is modeled using exponentials. The ratio of 
Z boson to non-resonant Drell-Yan production is fixed 
by the standard model. 

The (two) jets arising from Higgs boson decay should 
contain 6-flavored quarks (b jets), whereas background 
from Z-l-jets has relatively few events with b jets. To im- 
prove the signal-to-background ratio, two of the jets in 
the events from the Z-|-dijet sample are required to ex- 
hibit properties consistent with those of jets containing 
6 quarks. The same 6-jet identification algorithm is 
used for the dielectron and dimuon samples. It is based 
on the finite lifetime of b hadrons giving a low probabil- 
ity that these tracks appear to arise from the primary 
vertex and considers all central tracks associated with a 
jet. A small probability corresponds to jets with tracks 
with large impact parameter, as expected in b hadron 
decays. The efficiency for tagging a 6 jet from Higgs de- 
cay is approximately 50%, determined as described in the 
next paragraph. The probability of misidcntifying a jet 
arising from a charm quark as a 6 jet is roughly 20%. 
The probability to misidentify a jet arising from a light 
quark (u, d, s) or gluon as a & jet is roughly 4%. This 
choice of efficiency and purity optimizes the sensitivity 
of the analysis. The relatively large pcr-jet light-flavor 
misidentification rate can be accommodated because two 
tagged jets are required in each event. 

For background yields determined from simulated 
events, the probability as a function of jet pT and rj that 
a jet of a given flavor would be identified (tagged) as a 
b jet is applied to each jet in an event. The probability 
functions are derived from control data samples. For jets 
in the simulated events, the flavor is determined from a 
priori knowledge of the parton that gives rise to the jet. 
The probability of having two 6-tagged jets is defined by 
convoluting the per-jet probabilities assuming there are 
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FIG. 1: The dijet invariant mass distribution in double- 
tagged Z+dijet events. The Higgs signal corresponds to 
Mh = 115 GeV /(? . (The uncertainties are statistical only.) 



no jet-to-jet correlations introduced by the b-tag require- 
ment. The observed number oi Z + 2 6-jet events and the 
predicted background levels are shown in Tabic |T1 

The invariant mass of the two b jets in the Z + 2 b jet 
sample is shown in Fig. [1] This distribution is searched 
for an excess of events. The peak position in the dijet 
mass spectrum is expected to be at a lower value than 
the hypothetical Higgs mass because the jet energy is 
corrected to reflect the energy of particles in the jet cone 
without a general correction for the lower b jet response 
compared with light jets. If a muon is within AR < 
0.5 of the jet, then twice the muon momentum is added 
to the jet momentum. This is an approximation to the 
energy of both the muon and the accompanying neutrino. 
The expected contribution from Higgs boson production 
shown in Fig. [1] corresponds to Mh = 115 GeV/c^. 

Systematic uncertainties for signal and background 
arise from a variety of sources, including uncertainties on 
the trigger efficiency, on the corrections for differences 
between data and simulation for lepton reconstruction 
and identification efficiencies, lepton energy resolution, 
jet reconstruction efficiencies and energy determination, 
^-identification efficiency, uncertainties from theory and 
parton distribution functions for cross sections used for 
simulated events and uncertainties on the method used 
for instrumental background estimates. The uncertain- 
ties from these sources are shown in Table HH These are 
evaluated by varying each of the corrections by ±lcr, by 
comparing different methods (for the instrumental back- 
grounds), and by varying the parton distribution func- 
tions among the 20 error sets provided as part of the 
CTEQ6l library. The variations seen for different pro- 
cesses for a given uncertainty arise because of differences 
among the various background processes and because of 
intrinsic differences in the kinematic spectra from differ- 



ent Higgs mass hypotheses. 

The observed yield is consistent with background pre- 
dictions. Upper limits on the ZH production cross sec- 
tion are derived at 95% confidence level using the CLg 
method [1^, a modified frequentist procedure, with a log- 
likelihood ratio classifier. The shapes of dijet invariant- 
mass spectra of the signal and background are used to 
produce likelihoods that the data are consistent with the 
background-only hypothesis or with a background plus 
signal hypothesis. Systematic uncertainties are folded 
into the likelihoods via Gaussian distribution, with corre- 
lations maintained throughout. The data yield, predicted 
backgrounds and expected and observed limits are shown 
in Tablc Hnl for five hypothetical Higgs masses. The limits 
are also shown in Fig. [2l 

In summary, we have carried out a search for associated 
ZH production in events having two high-p^ electrons or 
muons and two jets identified as arising from b quarks. 
Consistency is found between data and background pre- 
dictions. A 95% confidence level upper limit on the Higgs 
boson cross section (7{pp ZH) x B{H bb) is set 
between 4.4 pb and 3.1 pb for Higgs bosons with mass 
between 105 GeV/c^ and 145 GcV/c^, respectively. 
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FIG. 2: The expected and observed cross-section limits are 
shown as a function of Higgs mass. The cross section based 
on the SM is shown for comparison. 
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TABLE I: Number of observed and expected background events. 







Z+>2 jets 






2 b tags 




Final state 


Z ee 


Z — > ^+/x~ 


Z — > £+£- 


Z — > cc 


Z — > jj,^ ^ 


Z — » £+£- 


Zbb 


9.1 


8.3 


17.4 


2.0 


1.3 


3.3 


Zjj 


414 


437 


851 


1.2 


2.6 


3.8 


tt 


2.7 


9.6 


12.3 


0.80 


3.1 


3.9 


zz + wz 


9.2 


21.4 


30.6 


0.32 


0.42 


0.74 


Instrumental 


28.0 


16.1 


44.1 


0.18 


0.41 


0.59 


Total background 


463 


493 


956 


4.5 


7.8 


12.3 


Observed events 


463 


545 


1008 


5 


10 


15 



TABLE II: Systematic uncertainty in background and signal 
predictions given as the fractional uncertainty on the event 
totals. The ranges correspond to variations introduced by 
different processes (background), the dijet mass window re- 
quirement (background and signal) and intrinsic differences in 
kinematics arising from different hypothesized Higgs masses 
(signal). 



Source 


Background 


Signal 


Lepton ID Efficiencies 


11% - 16% 


11% - 12% 


Lepton Resolution 


2% 


2% 


Jet ID Efficiency 


5% - 11% 


8% 


Jet Energy Reconstruction 


10% 


7% 


6-jet ID Efficiency 


10% - 12% 


9% 


Cross Sections 


6% - 19% 


7% 


Trigger Efficiency 


1% 


1% 


Instrumental Background 


2% (ee) 
12% (nfj.) 
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TABLE III: Numbers of predicted background and signal events and the observed yield after all selection requirements, including 
the addition of a dijet mass window. The mass window is centered on the mean of the reconstructed Higgs mass in simulated 
ZH events and has a width of ±1.5(t in which cr is the result of a gaussian fit to the reconstructed dijet mass distribution. The 
upper bounds differ slightly between the Z — > ee and Z events because of different resolutions. The window is applied 

for illustration, showing the yields in the region of highest predicted signal-to-background ratio. Also shown are the expected 
and observed upper limits on the cross section for the combined analysis at 95% confidence level computed as described in the 
text (without the mass window, but weighted by the bin-to- bin signal-to-background ratio). 





Mh ' 
ee 


= 105 GeV/c^ 
fifi 


Mh ■ 
ee 


= 115 GeV/c^ 
MM 


Mh ■ 
ee 


= 125 GeV/c^ 
fifi 


Mass window(GeV/c^) 


[65,113] 


[65,118] 


[72, 125] 




[70, 128] 


[75, 136] 


[78, 137] 


Predicted signal 


0.07 


0.06 


0.05 




0.05 


0.04 


0.03 


Background 


1.4 


3.1 


1.3 




3.1 


1.4 


2.8 


Data 


2 


3 


1 




3 


1 


4 


Expected (795 
Observed 0-95 




4.2 pb 
4.4 pb 




4.1 pb 
4.0 pb 






3.4 pb 
3.3 pb 





Mh 


= 135 GeV/c^ 


Mh 


= 145 GeV/c^ 




ee 


fifi 


ee 


MM 


Mass window(GeV/c'') 


[82, 143] 


[84, 147] 


[87, 156] 


[92, 160] 


Predicted signal 


0.027 


0.022 


0.015 


0.01 


Background 


1.6 


2.9 


1.6 


2.8 


Data 


1 


5 





6 


Expected 0-95 




2.8 pb 




2.6 pb 


Observed (J95 




3.1 pb 




3.4 pb 



